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I discuss some recent attempts to revive two old ideas regarding an 
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1. Introduction 



In this lecture I shall discuss some recent attempts to revive some 
old ideas to address the problem of solving QCD. I believe that it is 
timely to return to this problem which has been woefully neglected for 
the last decade. QCD is a permanent part of the theoretical landscape 
and eventually we will have to develop analytic tools for dealing with 
the theory in the infra-red. Lattice techniques are useful but they have 
not yet lived up to their promise. Even if one manages to derive the 
hadronic spectrum numerically, to an accuracy of 10% or even 1%, we 
will not be truly satisfied unless we have some analytic understanding 
of the results. Also, lattice Monte-Carlo methods can only be used to 
answer a small set of questions. Many issues of great conceptual and 
practical interest-in particular the calculation of scattering amplitudes, 
are thus far beyond lattice control. Any progress in controlling QCD in 
an explicit analytic fashion would be of great conceptual value. It would 
also be of great practical aid to experimentalists, who must use rather 
ad-hoc and primitive models of QCD scattering amplitudes to estimate 
the backgrounds to interesting new physics. 

I will discuss an attempt to derive a string representation of QCD 
and a revival of the large N approach to QCD. Both of these ideas have 
a long history, many theorist-years have been devoted to their pursuit- 
so far with little success. I believe that it is time to try again. In part 
this is because of the progress in the last few years in string theory. 
Our increased understanding of string theory should make the attempt 
to discover a stringy representation of QCD easier, and the methods 
explored in matrix models might be employed to study the large N limit 
of QCD. For both political and intellectual reasons I fervently urge string 
theorists to try their hand at these tasks. 
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2. QCD as a String Theory 



It is an old idea that QCD might be represented as a string the- 
ory. This notion dates back even before the development of QCD. In- 
deed, string theory itself was stumbled on in an attempt to guess sim- 
ple mathematical representations of strong interaction scattering ampli- 
tudes which embodied some of the features gleamed from the experi- 
ments of the 1960's. Many of the properties of hadrons are understand- 
able if we picture the hadrons as string-like flux tubes. This picture 
is consistent with linear confinement, with the remarkably linear Regge 
trajectories and with the approximate duality of hadronic scattering 
amplitudes. 

Within QCD itself there is internal, theoretical support for this idea. 
First, the expansion of weak coupling perturbation theory can be 
interpreted as corresponding to an expansion of an equivalent string 
theory in which the string coupling is given by . This is the famous 
result of 't-Hooft's analysis of the expansion of perturbative QCD 
[1]. The same is true for any matrix model-i.e. a model invariant under 
SU(N) or U(N), in which the basic dynamical variable is a matrix 
in the adjoint representation of the group. The Feynman graphs in 
such a theory can be represented as triangulations of a two dimensional 
surface. This is achieved by writing the gluon propagator as a double 
index line and tiling the graph with plaquettes that cover the closed 

index loops. 't-Hooft's principal result was that one can use to pick 
out the topology, i.e. the genus=number of handles, of the surface, since 
a diagram which corresponds to a genus G Riemann surface is weighted 
by (jj) 2G ~ 2 - The leading order in the expansion of the free energy in 

powers of is proportional to TV 2 (reasonable since there are TV 2 gluons, 
and is given by the planar graphs of the theory. 

Another bit of evidence comes from the strong coupling lattice for- 
mulation of the theory. The strong coupling expansion of the free energy 
can indeed be represented as a sum over surfaces [5]. Again there is a 
natural large N expansion which picks out definite topologies for these 
surfaces. This result is an existence proof for a string formulation of 
QCD. However, the weights of the surfaces are extremely complicated 
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and it is not at all clear how to take the continuum limit. 

From quite general considerations we expect that the large N limit 
of QCD is quite smooth, and should exhibit almost all of the qualitative 

features of theory. Thus an expansion in powers of | or (^) 2 might 
be quite good. The longstanding hope has been to find an equivalent 
(dual) description of QCD as some kind of string theory, which would 
be useful in to calculate properties of the theory in the infrared. 

The problems with this approach are many. First, if QCD is de- 
scribable as a string theory it is not as simple a theory as that employed 
for critical strings. It appears to be easier to guess the string theory 
of everything than to guess the string theory of QCD. Most likely the 
weights of the surfaces that one would have to sum over will depend 
on the extrinsic geometry of the surface and not only its intrinsic ge- 
ometry. We know very little about such string theories. Also there are 
reasons to believe that a string formulation would require many (per- 
haps an infinite) new degrees of freedom in addition to the coordinates 
of the string. Finally, there is the important conceptual problem-how 
do strings manage to look like particles at short distances. The one 
thing we know for sure about QCD is that at large momentum transfer 
hadronic scattering amplitudes have canonical powerlike behavior in the 
momenta, up to calculable logarithmic corrections. String scattering, on 
the other hand, is remarkably soft. Critical string scattering amplitudes 
have, for large momentum transfer, Gaussian fall-off [3]. How do QCD 

strings avoid this?^ 



★ There is also the problem that for large TV there is typically a phase transition 

between the strong and weak coupling regimes [2]. 
f Recently there have been some interesting speculations regarding this problem 

[4]- 



3. Two Dimensional QCD 



Two dimensional QCD (QCD 2 ) is the perfect testing ground for the 
idea that gauge theories can be equivalent to string theory. First, many 
features of the theory are stringier in two than in four dimensions. For 
example, linear confinement is a perturbative feature which is exact at 
all distances. Most important is that the theory is exactly solvable. 
This is essentially because in two dimensions gluons have no physical, 
propagating degrees of freedom, there being no transverse dimensions. 
In fact QCD 2 is the next best thing to a topological field theory. The 
correlation functions in this theory will depend, as we shall see, only on 
the topology of the manifold on which formulate the theory and on its 
area. For this reason we will be able to solve the theory very easily and 
explicitly. 

Consider for example the expectation value of the Wilson loop for 

pure QCD 2 , Tr R Pe$c A » dx " , for any contour, C, which does not inter- 
sect itself. Choose an axial gauge, say A\ = 0, then the Lagrangian is 
quadratic, given by ^Tri? 2 , where E = d\Ao is the electric field. The 
Wilson loop describes a pair of charged particles propagating in time. 
This source produces, in two dimensions, a constant electric field. The 
Wilson loop is then given by the exponential of the constant energy of 
the pair integrated over space and time. This yields, 

Tr^ Pe$c A ^ =e -9 2 C 2 (R)A ^ (31) 

where g is the gauge coupling, C 2 (l?) the quadratic Casimir operator for 
representation R and A the area enclosed by the loop. The expectation 
value of more complicated Wilson loops that do self intersect can also be 
calculated. Kazakov and Rostov worked out a set of rules for such loops 
in the large N limit [6] . They are quite complicated. QCD 2 with quarks 
is also soluble, at least in the large N limit. The meson spectrum was 
solved for N — > oo by 't Hooft. It consists of an infinite set of confined 
mesons with masses m n that increase as vr? n ~ n. This provides one 
with a quite realistic and very instructive model of quark confinement 

[1], [7]- 

Is QCD 2 describable as a string theory? The answer is not known, 
although there is much evidence that the answer is yes. I shall describe 
below a study that I have carried out to investigate this issue [8]. 
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To simplify matters I shall discard the quarks and consider the pure 
gauge theory. This would correspond to a theory of closed strings only, 
quarks are attached to the ends of open strings. We shall consider the 
partition function for a U (N) or SU (N) gauge theory, on an arbitrary 
Euclidean manifold Ai, 

Z M = J [VA"]e~& Im^^f^f^ ( 3-2 ) 

One might think that in the absence of quarks the theory is totally 
trivial, since in two dimensions there are no physical gluon degrees of 
freedom. This is almost true, however the free energy of the gluons will 
depend non-trivially on the manifold on which they live. In fact, one 
cannot simply gauge the gluons away. If, for example, contains a 

non-contractible loop C, then if Tr ^ 1, one can not gauge 

to zero along C. Thus, the partition function will be sensitive to the 
topology of A4. 

Although non-trivial the theory is extremely simple, almost as sim- 
ple as a topological theory. It is easy to see that the partition function 
will only depend on the topology and on the area of the manifold M.. 
This is because the theory is invariant under all area preserving diffeo- 
morphisms. To demonstrate this note that the two-dimensional field 
strength can be written as F^ u = e^ u f ', where is the anti-symmetric 

tensor and / a scalar field. Thus the action is S = / Tr f 2 dfi, where 
d/A = ^/gd 2 x is the volume form on the manifold. This action is inde- 
pendent of the metric, except insofar as it appears in the volume form. 
Therefore the theory is invariant under area preserving diffeomorphisms 
(Woo). The partition function can thus only depend on the topology 
and on the area of the manifold A4, 

Z M = Z[G,g 2 ,A,N} = Z[G,g 2 A,N} , (3.3) 

where G is the genus of M.. 

Now we can state the conjecture that the logarithm of this partition 
function, the free energy, is identical to the partition function on some 
string theory, with target space A4, where we would identify the string 
coupling with -k and the string tension with g 2 N ', 



7 



\n(Z[G, <? A, N] ) = 4^5 Space M fata = 1, a' = S 2 iV] . (3.4) 



As a candidate for the type of string theory I am thinking of consider 
the Nambu action, wherein 



Z S M ng = E (ftt^yiWKJe/^VS , (3.5) 

fo=genus 

where g is the determinant of the induced metric 

dx^ 1, Ox^ 

9 = det[g a(3 ] = det[——-—G flu (x)} , (3.6) 

and G ^11/ [x) is the metric on the manifold A4. This string theory, when 
the target space is two-dimensional, is indeed invariant under area pre- 
serving diffeomorphisms of the target space. To see this note that 

\fg = \^§z\VG, which is obviously unchanged by a map x^ — > 

as long as |J^| = 1* Unfortunately the only way we know to quan- 
tize this theory is to transform it into the Polyakov action, which upon 
quantization yields the standard non- critical string [10]. This is not 
what we want to do here, since the resulting theory is not even Lorentz 
invariant. Is there another quantization of the Nambu string that differs 
from the Polyakov quantization in two- dimensions? The answer is not 
known. 



* Actually the Polyakov action with a two-dimensional target space also has a 
Woo symmetry, although is is realized in a very nonlinear fashion. One might 
speculate that this is related to the well known symmetry of the c = 1 
string theory [9]. 



3.1. Evaluation of the Partition Function 



The partition function for QCD 2 can easily be evaluated by means 
of the following idea, originally due to Migdal [11]. The trick is to 
use a particular lattice regularization of the theory which is both exact 
and additive. For the lattice we take an arbitrary triangulation of the 
manifold and define the partition function as 

Z M = jY[dU L Y[Zp[U P ] , (3.7) 

L plaq. 

where Up = ElLeplaq Ul, and Zp[Up] is some appropriate lattice ac- 
tion. Any action will do as long as it reduces in the continuum limit 
to the usual continuum action. Instead of the Wilson action, Zp{U) = 

e a 2 Tr ( u+U ) ? we shall choose the heat kernel action, 

Zp = J2 d RXR(Up)e-° 2c ^ Ap , (3.8) 

R 

where the sum runs over representations R of SU(N) (or U(N)), dp is 
the dimension of R, Xr{Up) the character of Up in this representation, 
C2(R) the quadratic Casimir operator of R and Ap the area of the 
plaquette. 

It is easy to see, using the completeness of the characters to ex- 
pand about Zp Ul 1 tX A ' idx ^2 r oIrxr(Up) = 5{Up — 1) + . . ., that in 
the continuum limit of this theory reduces to ordinary Yang-Mills the- 
ory. What is special about the heat kernel action is that it is additive. 
Namely, we can integrate over each link on the triangulation, say U\, 
which appears in precisely two triangles, using the orthogonality of the 
characters, / dVxa{XV)xb{V ] Y) = ^Xa(XY), and obtain, 

J dU 1 Zp 1 (U 2 U 3 U 1 )Zp 2 (ulu 4 U 5 ) = Zp 1+ p 2 (U 2 U 3 U A U 5 ) . (3.9) 

This formula expresses the unitarity of the action, since in fact Zp(U) = 
(U\e~ 9 ^ A |l), where A is the Laplacian on the group. 
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We can use this remarkable property of the heat kernel action to 
argue that the lattice representation is exact and independent of the 
triangulation. This is because we can use (3.9) in reverse to add as 
many triangles as desired, thus going to the continuum limit. On the 
other hand we can use (3.9) to reduce the number of triangles to the bare 
minimum necessary to capture the topology of M.q- A two- dimensional 
manifold of genus G can be described by a 4G-gon with identified sides: 

a\bia^ 1 b^ 1 . . . clg^ , G c ^g'^c^ ■ The partition function can be written using 
this triangulation as, 



Zm g = Y. d R e ~ 92C2{R)A / W^UiVViXRp^uM . . . U G V G U G V G ] . 

R 

(3.10) 

We can now evaluate the partition function using the orthogonality of 
the characters and the relation, f VUxa[AU BU^} = ^-Xa[A]xa[B], to 
obtain [12], [13], 

Z MG = Y. d R 2G z- W) > (3-H) 
R 



where A = g N is kept fixed. Thus we have an explicit expression for 
the partition function. It depends, as expected, only on the genus and 
the area of the manifold. 



3.2. The Large N expansion 

The formula (3.11) for the partition function is quite complicated, 
being written as a sum over all representations of SU(N). The repre- 
sentations of SU (N) or U (N) are labeled by the Young diagrams, with 
m boxes of length n\ > n2 > n% > . . . n m > 0. Such a representation 
has, 
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m m 



1=1 i=l 

dR = -KW,' h ' = N + n '- 1 ' h° = N-i (3.12) 

aw= n ^~ h i) ■ 

\<i<j<N 

Thus we have a very explicit sum and one can, in principle, expand each 
term in powers of and evaluate the sum. 

What do we expect if the string conjecture is correct? Consider the 
expansion in powers of of the free energy, 
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\n[Z Ma ] = J^fg (AA) • (3.13) 

9=0 

If this were given by a sum over maps of a two-dimensional surface of 
genus g onto a two-dimensional surface of dimension G we would expect 

that f^(XA) ~ (^) 25 ' _2 e _A " 4n , where n is the winding number of the 
map, i.e. the topological index that tells us how many times the map 
x(£) covers A4. This is the integral of the Jacobian of the map £ — > x, 

J d 2 ^det[^r], which differs from the Nambu area, J d 2 ^\ det[^-]|, since 
the surface can fold over itself. 

Now there is a minimum value that G can take, given the genus G 
of the target space and the winding number n. Thus for example there 
are no smooth maps of a sphere onto a torus or a torus onto a genus two 
surface. Similarly there are no smooth maps of a genus g surface onto 
a genus g surface that wind around it more than once. To get an idea 
of the bound consider holomorphic maps, in which case the Riemann- 
Hurwitz theorem state states that 2(g — 1) = 2n(G — 1) + B, where B 
is the total branching number. In the case of smooth maps there seems 
to be the following bound [14], 

2(0-1) > 2n(G- 1) . (3.14) 



Thus if QCD 2 is described by a string theory we would expect that 
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fg( XA ) ~ { e - nXA u n JA) otherwise ^ ^ ' ( 3 ' 15 ) 



We can use these inequalities as tests of whether our conjecture is cor- 
rect. To do this we need to expand (3.11) in powers of jj . 

The hardest case is that of the sphere (G = 0), since the sum over 
representations blow up rapidly and it is not even evident that there ex- 
ists a tamed large N expansion. We can break up the sum in (3.11) into 
a sum over representations with n boxes in the Young tableaux since, 

for large JV, C{R n ) N ^°° N^ i n l = Nn. Thus, 



^o = E4 e -- C2(K) ^T£4e- nAi (l + .-) • (3-16) 

R fi R n 



To evaluate this we need to evaluate the following sum, YlfR n • This 

can be done using a method of discrete orthogonal polynomials [8] , yield- 
ing, 

Rn 

Then it follows that, 



2XAe~ 2XA 



Z G =0 - exp[-N 2 [ln(l - e ) + _ & _ XA)2 + . . . ] +0(N°) + . . .] . 

(3.18) 

Here there are no constraints implied by the inequality (3.15) , but the 
structure of the expansion is very interesting. 

The case of the torus, (G = 1), is some what simpler. One can easily 
derive that (for SU(N)) [8], 
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Z G=1 =Y,e~^ C ^ -> exp[-iV lnr / (-e- AA )+ 

i (3.19) 

n=l ab=n ab+cd=n 



where 77(0:) = n^li(l — • This is totally consistent with the bound 
g>l. 

Most interesting is the case of G > 1, where the inequalities are 
quite stringent. In this case one can easily derive [8], 

n r=rep of S n 

where the sum is over representations of the symmetric group S n and d r 
is the dimension of the r th representation of S n . Not only is this in total 
accord with our expectations, but one can also show that Ug(A) =, for 

g = 1 + n(G — 1), is precisely the number of topologically inequivalent 
maps on the genus g manifold onto the genus G manifold with winding 
number n[15]. 

So the large N expansion of QCD 2 looks precisely like what we would 
expect from string considerations. What remains to be understood are 
the all the rational numbers that appear as coefficients of the powers 

of e~ XA and of jj in terms of the counting of maps of M g onto M.g- 
Some of these are understood, but not all. Then it remains to construct 
a string action that reproduces these counting rules. 
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4. Induced QCD 

4.1. The Large N Limit of QCD 

QCD is hard to solve since it is a theory with no free, adjustable 
or small parameters. In pure QCD (no quarks) the only parameter we 
can adjust is the number of colors, N. Luckily, in the large N limit 
QCD simplifies enormously, and this limit remains the best hope to 
yield an exact or controllable treatment of the theory. We know that 
as N = oo only planar graphs survive. More generally we know that in 
terms of the appropriate variables the large N limit of gauge invariant 
observables is given, for N = oo by the master field, namely a solution 
of an appropriate classical equation of motion [16]. The large N limit is 
in the nature of a semi-classical expansion, with N playing the role of 
Planck's constant. Unlike the running coupling N does not vary with 
momentum and we expect the large N limit to be qualitatively correct 
for all momenta, to correctly capture the small distance asymptotic 
freedom of the theory as well as exhibit confinement at large distances. 
In the N = oo we should have an infinite spectrum of stable mesons 
and glueballs. Even baryons, bound states of N quarks, are describable, 
in this limit, as solitons of the effective Lagrangian for the master field 
[17]. Thus the hope has survived that we could find an exact solution of 
QCD for N = oo, which would yield the hadronic spectrum, and would 
be the starting point for a systematic large N expansion which could 
allow us to calculate scattering amplitudes. 

The standard method of solving a theory in the large N limit is 
to find an appropriate saddlepoint for the partition function. In the 
case of QCD this is difficult. Consider the standard (Wilson) lattice 
formulation of the theory, 

Z QCB = f]JVU L e-^^ Tr ^ UL+h ^ . (4.1) 

"* L 

The integrand behaves as the exponential of an action that is of order 
TV 2 , thus one might hope to evaluate it by saddlepoint techniques. How- 

A 7" 2 

ever, the measure is also of order c and therefore one must somehow 
get rid of N 2 degrees of integration before this can be done. The reason 
QCD is not yet solved in the large N limit is that no one knows how to 
reduce the theory to N variables per site. 
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Another theory which is also insoluble in the large N limit is the 
non-critical string with c > 1. Following the recent success of the matrix 
model solutions of string theory [18], we can construct such strings if 
we could deal with the large N limit of a scalar matrix model in D 
dimensions, say 

Z string = / > J^p^ e -NE l Trt/(^)+^E llA1= i...z 3 Tr(^ +At ) ^ 
J i 

which describes a scalar field on a D- dimensional lattice. The connec- 
tion with string theory is made in the usual way, the Feynman diagrams 
of the perturbative expansion of (4.2) correspond, in an expansion in 
powers of , to triangulations of two-dimensional surfaces. The scalar 
fields correspond to matter on this surface and thus, (4.2) , could yield, 
at the appropriate critical point where the mean number of triangles 
diverges, a c = D string theory. The standard approach to the large 
N- limit of such a theory is to diagonalize the matrices 0, i.e. to pass to 

radial coordinates, fa = f^A^J, where Xi is diagonal. In terms of these 
variables, 

^ tring =/n m * ra * A2 M 

(4.3) 

where Vij = fiifij and A(fa) = Yl i<:j (fa - <f>j). 

The next step is to integrate out the diagonalization matrices, Q^. 
We can change variables from the f^s, defined on the sites to the V^'s, 
defined on the links Ui^i = U<ij> ^ U^A* ~ Ul^j)- The 
constraints arise since the V^'s are pure gauge fields. If not for the 
constraints we could perform the integral over the Vifs and reduce the 
integral to one over TV variables per site that could be evaluated by 
saddlepoint techniques. It is these constraints that have prevented the 
construction of strings with c > 1. 

Now let us combine these two models to consider QCD with adjoint 
scalar matter, 
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•adj 

QCD - 




II VU L II V ^ G ~ 



L i 



(4.4) 



This theory is invariant under standard gauge transformations, 4>i - ► 

Vi4>iV^; U ^ — > ViUfj,V^ + , which allow us to diagonalize the 0's. However 

the presence of the Wilson action prevents us from handling this theory 
for large N. If set the gauge coupling to zero, we recover the previous 
model, since in this limit we can drop the Wilson action term, as long as 
we enforce the constraints, tr[£/p] = 1. However if we take the opposite 
limit, i.e. set g = oo, then we can simply drop the Wilson action and 
the model will be soluble in the large N limit. This is induced QCD 



Induced QCD has the one great advantage of being soluble, or at 
least reducible to a well defined master field equation. This is because 
the integral over the link matrices can now be performed. This is the 
famous Itzykson- Zuber integral [20] , 



This formula is very profound, underlies all the analysis of the c = 1, 
matrix model, and can be derived in many ways. One is the demonstra- 
tion that the integral is given exactly by the WKB approximation, and 
the answer is simply the sum over the N\ saddlepoints, for which are 
the U are permutation matrices. 

Although soluble this model appears to be very far from QCD, since 
asymptotic freedom instructs us to set the lattice coupling to zero, not 
infinity, in the continuum limit. However, Kazakov and Migdal argued 
that even though there is no kinetic term for the gauge field, it could 
be induced at large distances [19]. They argued that if one integrates 
out the scalar mesons (even in the case of noninteracting scalars with 
U($) = ^m 2 $ 2 ), then at distances large compared to a, one would 
induce in four dimensions an effective gauge interaction, 



[19]. 




(4.5) 
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Ses(U) ~ — -J H—22) ^ F lv + finite as a ^ . (4.6) 

This is simply the one loop vacuum graph for the scalars in a background 
gauge field, which is logarithmically divergent in four dimension. Now 

this looks very much like the ordinary Yang-Mills action, tr F^F^", 

if we recall that asymptotic freedom tell us that -^-^ = ln( M 2 Q2 ), 
where M g is a mass scale for QCD, say the glueball mass. We can there- 
fore identify these two (the fact that there are TV 2 scalars is crucial, as 
is the sign of the effective action which is due to the non-asymptotic 

freedom of the scalars.) If we do so then we find that, Mg = mil a 22. 
Thus in the continuum limit the adjoint scalars become infinitely mas- 
sive and decouple, but not before they have drive -\ up, from zero at 

distance a to the large QCD value atdistance ^, where >> ^ >> a. 

The basic idea is that the infrared slavery of the scalars, at the size of 
the lattice spacing, produces an effective gauge theory at a larger scale 
(much larger than the inverse scalar mass), which then produces the 
usual asymptotically free fixed point theory. 

There are many problems with this idea. For one the hard gluons 
are not absent and their contribution will overwhelm that of the scalars 
at short distances. Their asymptotic freedom is more powerful than 
the infrared slavery of scalars. Another issue is that the above theory 
possesses a much larger symmetry than the SU (N) gauge symmetry of 
the usual lattice action. It is not difficult to see that, in D dimensions, 
it is invariant under [D — 1) x [N — 1) extra local t/(l)-gauge symme- 
tries. This is because the transformation Ufj,(x)—> V / \(x)U fl (x)V^(x J rfia), 
leaves the action invariant as long as V^(x) is a unitary matrix that com- 
mutes with <$>(x). If Vn(x) were independent of \i then this would be the 
ordinary gauge invariance. Thus we have D — 1 new gauge symmetries, 
which are of course isomorphic to the special unitary transformations 
that commute with $ [21]. Thus V^{x) = D^(x)£l(x), where ft(x) is the 
unitary matrix that diagonalizes $ and D^(x) is diagonal. 

A subset of this symmetry is the, field independent, local Zn sym- 
metry, Un(x)—> Z^U^(x)zj l , where Z^ is an element of the center of the 
group. This symmetry alone prevents the Wilson loop from acquiring 
an expectation value. A Wilson loop contains different links, and thus 
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W(C) = (UlecUl) -> (Ul^)W(C) W(C) = 0. This symme- 
try must be broken if we are to recover the QCD fixed point from this 
formulation [22]. 

Finally, as we shall see, the simple Gaussian model is soluble and 
the answer is very simple and not equivalent to QCD [21]. However, 
there are interesting attempts to save the model and furthermore even 
if it does not yield a solution of QCD it might provide some interesting 
soluble matrix models which could yield new solutions of new string 
theories. Induced QCD is a matrix model and thus it corresponds to 
some kind of sum over surfaces. If we look at the Itzykson-Zuber integral 
we note that it could be expressed as, 

/(</>, X ) = exp IV 4? Tr x 2 + a IV </> 4 Tr X " + (TV 2 ) 2 (Tr X 2 ) 2 + . . .] . 

(4.7) 

These terms will affect the structure of the large N expansion of the 
Feynman diagrams, and can be interpreted as yielding extra weights 
when the two-dimensional surfaces intersect [23] . Thus this model corre- 
sponds, perhaps, to some kind of string theory with weights that depend 
on the extrinsic geometry. 

4.2. Solution of the Gaussian model 

To try to solve the model of induced QCD we first integrate out he 
Uls, then look for extrema of the effective action, 

Slfc] =^ 2 [^Tr^[/(^)+^^ln/(^,^)+^^lnA 2 (^ . 

i i,fi i 

(4.8) 

In the large N limit the integral will be dominated by a translation- 
ally invariant saddlepoint for the density of eigenvalues of the matrices 

p{x) = jj"}2a=iK x ~ 0a)- Migdal has derived the master field 
equation for the saddlepoint, using the Schwinger-Dyson equations that 
are satisfied by I((f),x) [24]. These are consequences of the fact that 
I satisfies tr[(-^ J^) fc ]i" = tr(x) k I- The net result is that one derives 

an equation for the function F(z) = f dzj^, whose imaginary part is 
lmF(u) = —7rp(v), 
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Re-F(A) = P / — ln[^ - = g^"* + ^]. (4.9) 

This equation is much more complicated than the usual Riemann-Hilbert 
problem that one obtains for simple matrix models. It is sufficiently 
non- linear and complex that one might imagine that it describes QCD. 

The master field equation simplifies dramatically for D = 1. This is 
because in one dimension the gauge field can be gauged away completely, 
thus the model is equivalent to a scalar field on a one-dimensional lattice. 
The large N limit of this model describes the c = 1 string on a discrete 
target space, a model which has been solved in the double- scaling limit 
for small lattice spacing [25] . It undergoes a phase transition at a finite 
lattice spacing and it might be very instructive to use (4.9) to explore 
this phenomenon. 

In particular for the quadratic potential the path integral is Gaus- 
sian, 

Z = l Y[V^ n e- N ^{^ 2 n-^n +1 } (41Q) 

J n 

Thus the eigenvalues of $ will be given by the semi-circular distribution, 

/ 2 2 

namely np(u) = \ I \i — where \i is determined by the mean of the 

squares of the eigenvalues, (-^Tr($ 2 )) = ^. It is therefore sufficient to 

calculate the expectation value of -^Tr($ 2 ), which is given by the one 
loop integral, 

I T r($ 2 ) = / = - 1 ee - . (4.11) 

N y J J 27rm 2 + 2cos/j Vm 4 - 4 V V ' 



N J 2tt m 2 + 2cos7J Vm 4 - 4 \i 

— TT 

It is easy to verify that this solves (4.9), using the fact that 



However, if we return to (3.8), we see that the integral involved is of the 
same form for any D, as long as $tV'{v) is linear in v. This suggests 
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that we can find a solution of (4.9) with a semi-circular distribution of 
eigenvalues for a quadratic potential in any dimension [21]. 

Indeed, one can see that a semi-circular distribution of eigenvalues 
satisfies (4.9) for any D as long as, 



, ±(D ) = " 2 (p-i)± 2 y_m«-4(2g-i) (4 13) 

This solution is much too trivial to describe QCD. In particular, for 
D > 1 there is no sensible continuum limit of the model. 

Is the solution unique? To see that it is note that in the master 
field equation the dimension of space-time enters only via the number of 

nearest neighbors of a given site, the coordination number of the lattice. 
The translationally invariant background scalar field is the same for any 
lattice with the same coordination number. The observables, say the 
scalar propagator, will of course depend on the full structure of the 
lattice, but not he background field. Therefore we can choose another 
simpler lattice with the same coordination number, say a Bethe lattice, 
which contains no closed loops. For such a lattice, as in the case of the 
D = 1 model, the gauge field can be eliminated completely, and the 
model is equivalent to, 

z^-fv*-"^**^™. (4 , 4) 

This model is easily soluble. We define Z((f>) to be the partition function 
of a branch of the Bethe lattice with coordination number 2D, so that 

Z = J D0Z(0) 2D e- I # iVTr ^ 2 . Z((f>) satisfies the equation, 

Z(<f>) = J V$Z^f D -\-^ N ^+ N ^^ . (4.15) 

These equations are easily soluble. Take Z((f>) to have the form Z{4>) = 

ce - N 2^ x $ ? then (4.15) determines a to equal a = m • 

Then (4.14) can be used to determine Tr 2 = m -2+ 2a D i wmcn agrees 
precisely with ^ as given by (4.13). 



* I thank C. Bachas for emphasizing this point to me. 
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4.3. Prospects 

The simplest Gaussian model fails, but all hope is not lost. It is 
certainly possible to induce QCD if one introduces enough flavors of 
matter. The problem is that one then loses solubility. It might be that 
the self interactions of the scalars could be adjusted to drive the theory 
towards the asymptotically free fixed point. This hope has been pursued 
with great vigor by Migdal, who has also considered adding fermions, 
not too many so that the model remains soluble, so as to break the 
Z/v symmetry [26]. Time will tell whether this will succeed. Even if 
it does not these model might yield a new class of interesting soluble 
matrix models which could teach us something about new classes of 
strings, perhaps strings that depend on extrinsic geometry. For this 
reason alone it is worth studying these models. 
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